Outflow and Infall in a Sample of Massive Star Forming Regions 

II: Large Scale Kinematics 



P. D. Klaassen, & C. D. Wilson 
Dept. of Physics and Astronomy, McMaster University, Hamilton, ON, Canada 

klaasspSphysics .mcmaster . ca 

ABSTRACT 



We present maps of seven sources selected from iKlaassen &: WilsonI (l2007al ) in 
SiO (J=8-7) and HCO+ and R^^CO+ ( J=4-3) which were obtained using HARP- 
B on the James Clerk Maxwell Telescope. We find that four out of our seven 
sources have infall signatures based on their HCO+ emission profiles. From our 
maps, we have determined the extent of both the outflowing and infalling regions 
towards these sources, and constrained the amount of infalling and outflowing 
mass as well as the mass infall rate for each massive star forming region. From 
our SiO observations, we estimate the source luminosity required to shock the 
surroundings of these massive star forming regions and find luminosities similar 
to those of the HII regions themselves. We find that the ratio between our infall 
and outflow masses is less than one, suggesting high mass entrainment rates in 
the molecular outflows. We also find that the large scale molecular infall rate 
towards GlO.6-0.4 is comparable to the small scale molecular infall rate found in 
previous studies. 

Subject headings: Stars: Formation - ISM: Jets and Outflows - HII regions - 
submillimeter - ISM: kinematics and dynamics 



Introduction 



In the last 20 years, the processes involved in the formation of solar type stars have 
been well characterized. An evolutionary sequence has been developed in which a prestellar 



Andre et al 



Andre et al 



20001) b egins to collapse under its own gravity and forms a Class 



19931 ) with a large disk of material surrounding it. This accreting 



core [I.e. 
object (i.e 

protostar begins to pro duce an observable molecular out flow which releases the build up of 
angular momentum (i.e. lBachilleiill996l : lArce et al.ll2007l ) and evolves into the Class 1 stage. 
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Later, as the protostar becomes more revealed, accretion slows and a pre- main- sequence star 
forms. It then undergoes quasi-static cont raction uritil co re hydrogen burning commences 
and the star reaches the main sequence (i.e. IShu et al.lll987l ). These processes have been well 



studied and characterized based in part on the high resolution observations possible towards 
the closest examples of low mass star forming regions. 



The initial mass function (IMF) of stars (i.e. lKroupall200ll ) shows that most of the mass 
contained in stars is within the lower mass stars, and that massive stars do not form nearly 
as often as their lower mass counterparts. There are a number of nearby low mass star 
forming regions (i.e. Taurus, Ophiucus, Chameleon, and Serpens are well studied nearby 
examples) in which there is no massive star formation. High mass sta rs do not, in genera l, 
form in isolation, but together with low mass stars (i.e. Orion, see iRiddick et al.l 120071 ). 
The high densities required for massive star formation also mean that more of them form 
in the molecular ring. This means that statistically, the closest examples of regions forming 
high mass stars are further away, and are also forming clusters of stars. Both distance and 
clustering can cause confusion when attempting to determine the dynamics in a region. 

Add to the large average distances and source confusion the higher energies involved in 
the formation of massive stars and the presence of HII regions surrounding the more evolved 
massive protostars, and it becomes clear why high mass star formation is not nearly as well 
understood as lower mass star formation. If we assume that massi ve stars form in a manner 
similar to their lower mass counterparts (i.e. iMcKee fc T^boosh . i t is possible to su ggest 



an evolutionary sequence starting from an Infrared Dark Cloud (i.e. iPillai et al.ll2006l ) with 



a deeply embedded protos tar which then begins to heat its surroundings and creates a hot 



core I.e. 



Kurtz et al.ll2000l ). The protostar then begins burni ng hydroge n brightly enough to 



towards and up the main sequence (IBeuther fc Shepherd 



2005 



ionize its surroundings, forming a hyperc ompact HII region (IKetd 120031) . which continues to 
expand into an ultracompact HII re gion (Wood &: Churchwell 1989) as the protostar evolves 



Keto fc Woodll2006f ). There 



are many observable phenomena, such as large scale infall and outflows, which appear to be 
common between both low and high mass star formation, which suggests that perhaps high 
mass stars may form via the accretion of matter onto a single object in a manner similar to 
low mass stars. However, for high mass star formation, there are still a number of questions 
to be answered. For the large scale phenomena observable with a single dish telescope, we 
still do not know how large an area the outflowing gas shocks, how much energy is required 
to shock these regions, how large the infall region is, or how much mass from the protostellar 
core undergoes infall onto the central massive object. 



Here, we present follow-up observations to iKlaassen fc Wilson! (j2007al . hereafter Paper 
1) with which we address the questions posed above for seven high mass star forming re- 
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gions mapped using the James Clerk Maxwell Telescop€0. In Paper 1, we obtained single 
pointing observations towards 23 raassive star forming regions previously known to have HII 
regi ons (Wood fc Churchwell 1989 : Kurtz et al.lll994jl and mo lecular outflows as traced by 



CS Jpiume et aljll997l : lshirlev et al 



20031 ) or CO flHuntei]ll997h . We were able to show that, 



of our observed high mass star forming regions with ongoing outflow activity (14 sources, as 
traced by SiO emission), half (seven) showed evidence for large scale infall motion as well 
(as traced by HCO"*" line asymmetries). We suggested that the sources with infall and out- 

1 actively accreting material (possibly in an ionized form, like that 



flow detections were sti 



discussed in iKetd (120071). while the other half may have already flnished accreting (like the 



source discussed in iKlaassen et al.ll2006l ). However, with single pointing observations, we 
could not determine the sizes of the infalling or shocked (outflowing) regions, and therefore 
could not discuss properties such as the extents of the infall and outflow regions, whether 
these phenomena were beam diluted, or how their sizes compared. In this paper, we present 
maps of seven sources in the same tracers as Paper 1 in order to understand better the nature 
of the large scale infall and outflow towards massive star forming regions. 

For this paper, we chose to map the SiO (J=8-7) and HCO^/H^^CO^ (J=4-3) emission 
towards seven massive star forming regions. Six sources were chose from Paper 1, and were 
chosen to be representative of the survey as a whole. Two sources (G29.96-0.02 and G5.89- 
0.39) showed only outflow signatures, while two showed both infall and outflow signatures 
(GlO.47+0.03 and G19.61-0.23), and one source (K3-50A) showed only a weak outflow sig- 
nature. It was unclear whether the emission towards the sixth source (G20.08-0.14) was due 
to one cloud (as suggested by the single peaked SiO line proflle) or two (as suggested by the 
double peaked HCO^ and H^^CO^ line proflles). The G20. 08 region was previou sly known 
to have two ultracompact HII regions in close proximity (i.e. ISewilo et al.ll2004l ). although 
the two HII regions are separated by more than a single JCMT beam. To this list of six 
sources, we added one more source, GlO.6-0.4. This source was added to determine whether 



the ionized and molecula r accretion observed on small scales in this region (i.e. lKeto fc Wood 
20061 : ISoUins et al.lboosh has a corresponding larger scale molecular infall signature. G10.6- 
0.4 will act as a good test for determining whether our large scale results reflect the smaller 
scale dynamics, since the small scale structures are fairly well understood for this region. 
From this point forward, all source names based on Galactic coordinates will be truncated 
to latitudes only. 



In Section [2] we present our observations and in Section [3] we determine outflow and 



^The James Clerk Maxwell Telescope is operated by The Joint Astronomy Centre on behalf of the 
Science and Technology Facilities Council of the United Kingdom, the Netherlands Organisation for Scientific 
Research, and the National Research Council of Canada. 
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infall properties based on these observations. Our results are then discussed in Section HI 
and we summarize in Section [5] 



2. Observations 

Observations of SiO (J=8-7), HCO+ and Hi3C0+ (J=4-3) were obtained at the JCMT 
towards G10.47, G10.62, G19.61, G20.08, G29.96 and K3-50A in April and August 2007 as 
project M06BC08 along with HCO+ observations of G5.89. SiO (J=8-7) an d H^^CO+ were 



obser ved simultaneously towards G5. 89-0. 39 as part of project M03AC13 (see iKlaassen et al. 



20061 . for a description of these observations) and taken from the JCMT archive for this study. 
The archive map is smaller than our new maps; however, both the SiO and H^^CO^ emission 
falls below the rms limits within the edges of this map. 

The spectral resolution of the new datasets is 0.42 km s~^, the observations are spaced 
at 7" to Nyquist sample the 14" full width of the JCMT beam at these frequencies, and the 
main beam efficiency was 0.62. Map center positions and local standard of rest velocities 
are given for our six original sources in Paper 1. For G10.6, these properties are as follows: 
a = 18:10:27.8, 6 = -19:56:04, and Vlsr = 3 k m s~^. For reason s presented in Section WT\ 



in this paper we adopt a distance of 5.8 kpc fi.e.lGibb et al.ll2004r) for G10.37, which is much 



closer than the value of 12 kpc (i.e. IShirley et al.ll2003l ) used in Paper 1 



The 2007 observations were obtained using HARP-B and the ACSIS autororrelator 
system in either grid position switching or jiggle position switching mode. Both methods 
provided similar noise limits in similar times. The noise levels in each map are given in 
the respective figure captions and were generally consistent across each map, staying within 
approximately 30% of a mean value when receivers with bad noise levels were removed from 
the maps. These observations were reduced using the Joint Astronomy Center Starlink 
software. The 2003 observations were obtained with RxB in raster mapping mode using 
the DAS autocorrelator system and reduced using SPECX. First order linear baselines were 
removed from each pixel within each map, and the reduced maps were exported to the 
MIRIAD format for analysis. 



3. Results 

SiO and HCO"*" were chosen for this study because they trace outflow and infall, respec- 
tively, in star forming regions on all mass scales. Once each map was properly calibrated and 
exported into MIRIAD, the peak brightness temperature and line full width at zero power 
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(FWZP) were determined from the spectrum at the map pointing center. The uncertainty in 
the integrated intensity for each hne was determined using A/ =Trms ^'^(Nchan)^^^ where Trms 
is the single channel rms uncertainty, Az/ is the velocity resolution of our observations, and 
Nchan is the number of channels used to cover the full width of the line (FWZP). Note that 
for sources with self absorbed HCO+ profiles, the peak line temperature is an underestimate 
of the source brightness temperature. 

In our maps, we find that the HCO^ emission towards each source is extended, as is 
the H^^CO"'" emission; however, H^^CO^ does not reach the same spatial scales as the 77 
times more abundant (Wilson & Rood, 1994) HCO"^. In only five of our seven sources is 
the SiO emission at least marginally resolved. The deconvolved source size was determined 
using Dsou = ^/ -Dobs ~ -^bcam whcrc the longest axis of the observed region full width at half 
maximum diameter (-Dobs) was used to compare to the beam diameter (-Dbeam)- Towards 
the two unresolved sources (K3-50A and G19.61), the deconvolved diameter of the emitting 
region is less than 50% of the JCMT primary beam, while the deconvolved diameter is greater 
than 70% of the JCMT primary beam towards our other sources. Our SiO detection is weak 
in K3-50A (the line peak was only detected at 3.5-4a) and was insufficient for detecting SiO 
beyond the central pixel. This was not the case for G19.61 (the line peak was detected at 
7a), yet the SiO emission in this source was also unresolved. Thus, if the shocked region 
in G19.61 were as extended as our JCMT beam, we would have detected the extended SiO 
emission. 

Emission maps for each source are shown in Figures [IE In each case, IICO+ emission 
spectra are shown for pixels with HCO"'" integrated intensities which are lOcr or greater. For 
each tracer, the intensity scale for the HCO"'" spectra is constant for a given source, and the 
peak HCO"*" intensity for each source is listed in Table [2l The dashed contour in each plot 
represents the 3a emission region for H^'^CO^ and the solid contours represent SiO emission, 
starting at 3a and increasing in intervals of 3a. The 3a rms uncertainties in the integrated 
intensities are given in each figure caption. The star in each figure marks the position of the 
peak in the HCO"*" integrated intensity. In all cases except K3-50A and G29.96, the positions 
of the HCO"*" peak is within the same beam as the SiO peak. For the two remaining sources, 
this offset may be due to pointing errors between the two maps. 

Using the FWZP from the central pixel, the integrated intensity for each species across 
each map was determined. The rms uncertainty in the integrated intensity was determined 
using the same method as presented in Paper 1. For each tracer, the line is broadest in the 
map center, and so using the same FWZP across the entire map captures all of the emission 
from that transition. The integrated intensities averaged over each map in SiO and HCO"*" 
are shown in Tables [1] and [21 respectively. 



- 6 - 



The integrated intensities in SiO and HCO^ were used to find the column density of 
gas emitting across the map in that tracer. For optically thin gas, the column density in 
the higher energy level is proportional to the integrated intensity of the line, with a scaling 
constant proportional to the frequency of t he transition a nd the Einstein coefficient for that 
transition (see Equation 1 of Paper 1, or iTielend l2005l ). Assuming local thermodynamic 
equilibrium and an ambient temperature of 44 K (as was done in Paper 1), we can use the 
partition function to determin e the column density o f optically thin gas in each tracer (see 
for instance Equation 14.38 of lRholfs fc Wilsonll2004l ). The average column densities of SiO 
and HCO^ are presented in Tables [1] and [21 The values for the H^^CO"*" column density are 
not presented here, but were found to be only 4-22 times less than those for HCO"'". That 
we do not see a larger difference in the column densities, given the abundance ratio of 77 
between the isotopologues, suggests that we are not seeing all of the HCO"*" emission because 
the line is optically thick. Thus, our derived column densities are lower limits to the true 
HCO"*" column densities due to optical depth effects. This trend was also seen in Paper 1, 
and supports the claim that the double peaked HCO+ profiles towards four of our sources 
are due to self absorption of the optically thick line, and not multiple clouds along the line 
of sight. 



3.1. Outflow properties derived from SiO observations 



In an environment which has recently been shocked, the abundance of SiO can be en- 
hanced up to a factor of 10~* (from dark cloud abundances of approximately 10^^^) with 



Caselli et al. 


1997; 


Schilke et al. 


1997) 



bining with oxygen to form SiO (i.e. 
imately 10^ years, SiO deplet es out of the gas phase again and the abundance drops (i.e. 



Pineau des Forets et al.l 119971 ). Because of the short enhancement lifetime, SiO is generally 
only detectable in regions which have recently been shocked. Where we detect SiO (J=4-3) 
we suggest there is ongoing outfiow activity. 

We present the column density of SiO towards each source in Table [1] We can divide 
this value by the abundance of SiO with respect to H2 in order to determine the column 
density of shocked gas in each region. These values are also presented in Table [H We used 
the same methods as described in Paper 1 to calculate the abundance of SiO towards G10.6 
(log[Ai'gio] + log[Xcs] - log[Ai'cs] = log[Xsio] = -9.92) and we have assumed a constant SiO 
abundance across each map. We note that the SiO column densities listed in Paper 1 are 
over estimated by a factor of 3.8, and thus the logarithmic SiO abundances listed in that 
paper are corrected in Table [H Because we have adopted a new distance for G10.47, we 
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have had in addition to modify our origina l abundance estimate from Paper 1. The CS 
column density quoted in lPlume et al.l (119971) appea rs to be distance independent; however, 
the CS abundance presented in IShirley et al.l (120031 ) is dependent on the assumed distance 
to the source (see their Equations 7 and 12). Thus, our SiO abundance has changed from 
-9.53 in Paper 1 to -10.42 when the new distance and revised column density are taken into 
account. The SiO abundance determined for K3-50A in Paper 1 (log[Xsio] = -11.6) is much 
closer to that found in dark clouds than in outflow regions. The region in which SiO is 
emitting is unresolved by our observations and highly beam diluted. Thus, for our analysis, 
we have assumed an SiO abundance of -10.15, the average abundance from our other six 
sources. Based on the abundance estimates presented in Table [H we have determined the 
total column of shocked material towards each source, from which we determined the mass 
of shocked gas because we know the area of the emitting region (see below). For our sources 
which are unresolved in SiO, we have taken the radius of our primary beam to be the radius 
of the emitting region. 

Comparing the integrated intensities of our SiO emission in the central pixel of each 
map, where our signal to noise ratio is highest, to the models presented in ISchilke et al. 



(Il997l ). we suggest shock velocities of 20-22 km s \ and ambient d ensities of 



10^ cm-3. 



The ambient densities suggested here are consistent with those of iPlume et al.l (119971 ) who 
determined densities and CS column densities for these regions using large velocity gradient 
models. For each of our sources, we adopt a shock velocity of Vs = 21 km s~^. 

Because we have mapped the SiO emission towards seven sources, we can determine the 
radius of the region in which SiO is being emitted (Rs). The method used to determine the 
deconvolved source size was described in Section [21 and the radius of the emitting region 
is presented in Table [TJ From this radius and the shock velocity determined above, we can 
constrain the ages of our outflows from simple kinematic arguments. These approximate 
outflow ages are presented in Table [2] and are used later to constrain the amount of infalling 
material. 

In order for the outflow shocks to reach their observed sizes given the mass of shocked 
gas, the powering sources for these outflows must produce enough energy to push through the 
ambient medium. If we assume a stea d y win d is driving a spherical bubble of gas outward, 
we can use the model of ICastor et al.l (119751) to de termine the luminosity required to push 
the ambient gas. Equation 6 of ICastor et al.l (119751 ) can be re-expressed as: 



^2 ^/1547r\ 



where Lq = Eot is the luminosity imparted by the wind, Rg is the radius of the bubble, Vg is 
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the shock velocity, and related to the outflow age by Vg = Rs/t, and po is the ambient density, 
and related to the ambient number density by the mass of molecular hydrogen [po = pmnno)- 
While much of this wind luminosity goes into accelerating the ambient medium, some of it 
is also transferred into thermal energy. From lBasu et al.l (119991 ). the thermal energy can be 
expressed as Eth = (5/ll)Lot, or, given the value for t above, Eth = {5/ll)LoRs/vs- Both 
the calculated wind luminosity and thermal energy required to power these wind driven 
shocks are listed in Table [H These thermal energies are of approximately the same order of 
magnitude as the kinetic energies calculated for the CO outflows in a number of massive star 
forming regions (see for instance IWu et al.l l2004l : iKlaassen &: Wilson! l2007bl ): similarly the 
wind luminosities are cora parable to the luminosities of their HIT regions (see for instance 



Wood fc Churchwelllll989[ ). 



3.2. Infall properties derived from HCO^ observations 

Many authors suggest that an asymmetric self absorbed line profile of an optically 
thick line tracer in which the blue emission peak is brighter than the red, and the optically 



thin tracer has a single line peak, can be indicative of infall (i.e. iMardones et al.l Il997 



Gregersen et al.l 119971 : iFuUer et al.ll2005l ). Here, our optically thick line is HCO+ (J=4-3) 



and the optically thin isotopologue is H^^CO"^ (J=4-3). 

From Paper 1, we expected to observe infall signatures in G10.47 and G19.61. Infall is 
detected towards these two sources, as well as towards G20.08 from our original study. In our 
single pointing observations, G20.08 was shown to have a double peaked line profile in both 
HCO"'" and H^'^CO^, which was suggested to be due to two overlapping sources. However, 
in the more extended map, we see that H^^CO"*" becomes single peaked, suggesting that 
perhaps it is optically thick towards the center of the source and that mapping is required 
to resolve the ambiguity with respect to the chance super-position of line of sight clouds 
and self absorption of multiple isotopologues. G10.6 was not in our original sample, but was 
added because it was known to have small scale infall motions. It too is seen to have a large 
scale infall signature. 

In Paper 1, we suggested that beam dilution may have been a contributing factor as to 
why we did not detect infall in a number of sources. That we have detected extended infall 
in two distant (~ 6 kpc) sources in this study suggests that if beam dilution of the signal 
is a problem, it only affects intrinsically smaller sources, and not intrinsically more distant 
ones, since our unresolved detections are for G19.61 and G20.08 at 4.5 and 4.1 kpc. 

Infall velocities were calculated for each self absorbed HCO+ spectrum towards G10.47, 
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G10.6, G19.61, and G20.08 using the models of iMyers et al.l (119961 ) in the same way as 
described in Paper 1. From our rms uncertainties in hne temperatures and the spectral 
resolution of our observations, we determined the uncertainties associated with our derived 
infall velocities. We used these uncertainties to limit the regions in which we claim infall 
detections. Only if the derived infall velocity is more than three times our uncertainty for 
that spectrum do we claim to detect infall in that pixel. Towards G19.61 and G20.08, robust 
detections were made in one and three pixels (respectively) suggesting that the infall in both 
of these sources is unresolved. Towards G10.47 and G10.6 (both at approximately 6 kpc), 
robust infall detections were made in seven Nyquist sampled pixels per map, suggesting that 
the large scale infall in these regions is extended. 

Along the northern, blue shifted, outflow lobe in G10.47, the HCO+ line profiles again 
become double peaked. It is quite possible that there is a second instance of a collapsing 
region here since SiO is marginally detected (at the 3a level) in the middle of this northern 
extension. The double peaked profiles do lead to infall velocities (in three pixels) which are 
greater than 3a; however, we do not include this region in our infall analysis for a number 
of reasons. The infall signature is not continuous from the central region to this northern 
region, which suggests that there are two distinct regions of infall. This discontinuity would 
confuse conclusions about the central infall region. Similarly, we are interested in the infall 
activity surrounding HII regions, whe re we know massive s tars are forming. Neither the mi d 
infrared nor radio maps of this region (jPascucci et al.l (120041 ) and lWood fc Churchwelll (jl989l ). 
respectively) extend far enough to the north to determine whether there are HII regions in 
this northern ext ension. Archival SCUBA maps of this region at the same resolution as 
our observations (jPi Francesco et al.l l2008l ) show emission which has the same extent and 
morphology as our HCO^ observations. 

For each region in which infall was observed, we determined an average infall velocity 
by taking the mean of the infall velocities determined for each position in each map using 
the models of iMyers et al.l (119961 ). This averaging resulted in the velocities shown in Table 
[2j We chose to use this method over creating an averaged spectrum and determining infall 
velocities from it. The latter method resulted in unreasonably high infall rates (i.e. 5 km 
s~^ in G10.47), which give very large mass infall rates (0.22 Mq yr~^) and total infall masses 
(~ 2000 Mq) assuming constant infall rates over the lifetime of the outflow. Averaging the 
spectra and then determining the infall velocity increased the infall velocity by increasing 
the full width at half maximum (FWH M) of the line. Sin ce the infall velocity varies as 
the FWHM squared (see Equation 9 of iMyers et al.l Il996l ). this broadening leads to an 
overestimate of the infall velocity. Thus, for our unresolved sources, the calculated infall 
velocities may be upper limits, since there could be unresolved velocity gradients within the 
JCMT beam. Higher resolution imaging of these sources could decrease our calculated infall 
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velocities. 

Using the derived average infall velocities for each region, we determined mass infall 
rates for each sources using M = dM/dt ^ pVvi^i/r = (4/3)7rno/imHr^t>in, where p = pmnno 
is the mass density, fin is the infall velocity presented in Table [2], and V/r is the area of the 
infalling area (which we assume to be circular). Since the area in which we detected an infall 
signature was always similar to the area in which SiO was detected, we used the radius of 
the SiO emitting region as the radius of the infalling region. 

From the radii of the SiO emitting regions and a shock velocity of 21 km s~^, we were able 
to estimate outflow ages in Section 13. 1[ Assuming infall and outflow are contemporaneous, 
these outflow ages can be used to constrain how much mass has undergone infall. These 
values are also given in Table [21 and represent the total mass of infalling molecular gas, not 
just the mass of the infalling HCO"'". 



Discussion 



It is interesting to note that the infall and outflow signatures in G19.61 are both unre- 
solved, and the signatures in G20.08 are both marginally resolved. This comparison suggests 
that the shocked region, as traced by SiO, may have a similar physical size as the infall re- 
gions in these smaller sources. Towards G10.6 and G10.47, the regions over which we have 
robustly detected extended infall, our deconvolved SiO emitting regions also appear to be 
on the same approximate spatial scales. 

Comparing infall mass and outflow mass, we note that our infall masse s are lower than 
our outflow masses. This is contrary to the results of iBeuther et al.l (l2002n : however, their 
ratio was calc ulated for the protos tellar jet, not the entrained molecular outflow. This 
comparison to iBeuther et al.l (|2002| ) suggests that more of the ambient mass is entrained 
into the large scal e molecular outflow (produced by a protostellar jet) than falls in onto the 



forming stars (i.e. IShepherd fc Churchwelllll996l ) 



4.1. G10.47: Extended outflow or chance super-position? 



The large spatial extent of the HCO+ emission in our map of G10.47 places it among the 
largest molecular outflows due to star formation in our Galaxy. Some Herbig-Haro objects 
have been observed emanating from nearb y lower mass star formin g regions at distances of 
up to 12 pc from their powering stars (i.e. iBally fc ReipurthI |2002| ) : however in those cases, 
the associated molecular outflow is not nearly as wide as the case of G10.47, and those jets 
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are not pushing through an ambient medium dense enough for massive stars to form (i.e. 
10^ cm"^). In Paper 1, we used a distance estimate of 12 kpc to G10.47. Given the extent 
of HCO"^ gas (~ 2') in our current maps, this distance would give a hnear scale of over 4 
pc f or a single outflo w. Other authors have suggested a distance to this source of 5.8 kpc 



I.e. 



Gibb et al.ll2004j ). For the analysis presented here, we adopt this new distance, which 
changes the length of the outflow from 4 pc to 2 pc. Given that this source is along a line of 
sight reasonably close to both the Galactic plane and Galactic center, it is not unreasonable 
to suggest that what appears to be an extremely large outflow is instead the super-position 
of a number of clouds along the line of sight. However, based on the arguments presented 
below, we suggest that the gas kinematics and morphology are due to a large outflow. 

Examining the spectra of our HCO+ map, we note that while the HCO+ emission is 
extended, it is not present in all pixels, suggesting that we have captured all of the emission 
in this velocity range. Figure [8] shows the HCO^ line profiles averaged over the blue and red 
shifted emission regions as well as the HCO^ and H^^CO"*" line profiles averaged over the 
central nine pixels of our map. The red and blue spectra, while having very little overlap 
between their profiles, do both fit within the width of the source center line profile. The 
velocities at line peak for the blue and red spectra are -64 and -70 km s~^ respectively, and 
match the velocities of the two peaks of the self absorbed, averaged HCO"*" central profile. 
The single peak of the optically thin H^^CO"'" is at -67 km s~^ as is the minimum of the 
HCO"*" self absorption feature. This self-absorption lies directly between the peaks of the red 
and blue profiles, and is the same as the source rest velocity. No other sources in the W31 
region have rest velocities similar to the ones seen in the G10.47 region. 

The double peak in the HCO+ spectrum is due to self absorption, not multiple compo- 
nents, since there is a single peak in the H^^CO"*" spectr um. The same is true for the CO 



isotopologues in this region. The BIMA observation s of iGibb et al.l (120041 ) show one peak 



in C^^O, while the ^^CO observations of lOlmi et al.l (119961 ) with the IRAM Interferometer 
show two peaks indicative of self absorption. 

The first moment map of the HCO"'" emission (grayscale in Figure [9]) in this region 
also suggests a single outflow. The observed smooth transition from blue to red along the 
source would be more likely to be observed in one source with a velocity gradient than three 
separate sources with discontinuous cloud boundaries. The second moment map of this region 
(contours in Figure [9]) suggest that the largest velocity dispersion is towards the center, and 
continuously decreases towards the emission edges (outflow lobe ends), suggesting that most 
of the mass is centrally concentrated. 



While we favor the outflow interpretation of this morphology, higher resolution, large 
scale observations of this source are necessary in order to resolve this ambiguity. This 
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region contains four individual UCHII regions and what is likely a large scale outflow with 
a coUimation factor (-RmajZ-Rmin) of approximately 1.5, which may be even more highly 
coUimated at higher resolution. Objects like G10.47, with large scale, fairly coUimated 
outflow emanating from a cluster of HII regions pose interesting challenges for accretion 
models. 



4.2. Large and Small Scale Infall in G10.6 



As stated previously, G10.6 was added to our source list because it was previ ously known 
to s how evidence for sma ll scale infall, both in molecular gas ( jSoUins et al.ll2005l ) and ionized 
gas (IKeto fc Woodll2006l ). These small scale infall signatures appear to be consistent with a 
single accretion flow which changes from being molecular to ionized as it passes the ionization 
boundary of the HII region. While we cannot determine the orientation of the large scale 
infall from our observations, we note that our derived mass infall rate (0.03 Mq yr~^) is 



comparable to, but slightly larger than, that of ISoUins et al.l (120051 ) (0.02 Mq yr^^). This 
similarity may suggest that the large scale infall possibly maps to the small scale accretion 
flow; in addition, despite resolving the infall region, we may still be over estimating the infall 
velocity due to velocity gradients in our larger JCMT beam (see Section ??). 



4.3. How mapping has changed the results from Paper 1 

None of our basic conclusions drawn from our single pointing observations of SiO (i.e. 
column densities in the central pixel or outflow lifetimes) have changed due to mapping save 
for the infall detection in G20.08. However, with respect to the outflowing gas, we have been 
able to place constraints on the energy required to shock these regions. With respect to 
infalling gas, we have placed better constraints on infall velocities and shown that the lack 
of infall signatures in some sources is not due to the signature being beam diluted out of 
observability at large distances. 

We have found that one of our sources which was previously not thought to have an 
infall signature (G20.08) in fact does appear to have large scale infall motions. This discovery 
suggests that more accurate infall results are possible with maps of massive star forming 
environments than with single pointing observations. Returning to the results of Paper 1, 
we determined that seven out of the 14 sources with ongoing outflow signatures had infall 
signatures as well. These numbers can now be revised to eight out of 15 outflow sources with 
infall signatures when G10.6 is added to the source list. Thus, we flnd that approximately 
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half of the sources with ongoing outflows have infall signatures which we take to trace 
accretion. Because our sources were selected based on the presence of HII regions, ongoing 
infall suggests that there is a phase of massive star formation in which the star is still actively 
accreting, but is also bur ning hydrogen brightly enough to ionize its surroundings (see for 
instance the discussion in lKeto fc KlaassenI 120081 ) . 



5. Conclusions 

We have found that just over half of the sources in our sample have large scale infall 
motions suggesting that the protostars powering the UCHII regions and massive outflows 
are still actively accreting material. Where we have detected infall motions, they have been 
on quite large scales, suggesting that infall, as traced by HCO+, is in general observable with 
single dish maps at the distances to massive star forming regions. Maps are required in order 
to capture all of the infalling gas since the infall signature can be extended, even at 6 kpc. 
In mapping these regions, we have been able to determine how much mass is falling into 
the central regions of our sample sources, which we would not have been able to determine 
with our single pointing observations. We find that when the infall is resolved, the infalling 
mass is greater than 30 Mp). Our derived m ass infall rate towards G10.6 is similar to the 



small scale rate found by lSoUins et al.l (120051 ) . suggesting that we may be observing the same 



physical process, but on different spatial scales. 

In general, we have found that SiO is not very extended; however, since we mapped 
these regions in SiO, we were able to determine the radius of the emitting region, the mass 
of shocked gas and how much energy is required to shock the gas surrounding a massive 
protostar. We find these energies and luminosities to be similar to the kinetic energy imparted 
to the outflowing gas and the luminosities of the HII regions themselves. We also find that 
the region over which we detect an infall signature and the region in which we detect SiO 
is quite similar, and find that in these regions, the total infall mass is less than the total 
outflow mass (Min/Mout ~ 0.1 to 0.8) assuming the two processes are contemporaneous. 



We would like to acknowledge the support of the National Science and Engineering 
Research Council of Canada (NSERC). P. D. K would also like to thank J. Wadsley for 
helpful discussions during the preparation of this manuscript. 



Table 1. SiO Observations and Derived Outflow Properties 



Source 


Dist 


Tmb 


FWZP^ 


X[SiO] 


<lT^y>dv> 


<log[Nsio]> 


<log[NHj> 


Mass 


Radius'^ 


L 


Eth 




(kpc) 


(K) 


(km s~^) 




(K km s-^) 


(cm~^) 


(cm~^) 


(Mo) 


(xlO^'^ cm) 


(xl03 Lo) 


(xlO^« erg) 


G5.89-0.39 


2.0 


2.37 


105 


-9.48 


4.38 


12.81 


22.29 


3.7 


1.9 


1.1 


0.18 


GlO.47+0.03 


5.8 


1.24 


30 


-10.42 


2.70 


12.60 


23.02 


184 


5.8 


10.3 


5.1 


GlO.6-0.4 


6.0 


1.14 


20 


-9.92 


8.74 


13.11 


23.03 


497 


9.4 


27.2 


21.6 


G19.61-0.23 


4.5 


0.71 


45 


-10.32 


4.57 


12.83 


23.15 


162 


<4.7 


<6.8^ 


<2.7 


G20.08-0.14 


4.1 


0.66 


30 


-11.06 


2.68 


12.60 


23.66 


374 


4.0 


4.8 1 


1.6 


G29.96-0.02 


9.0 


0.63 


40 


-9.70 


5.14 


12.88 


22.58 


124 


7.9 


19.2 


12.8 


K3-50A^^ 


8.6 


0.18 


20 


-10.15 


3.15 


12.67 


22.82 


276 


<9.0 


<21.7 


<18.7 



^Full Width at Zero Power of the SiO line. 

^Longest axis of SiO emitting region, after deconvolution from the 14" JCMT beam. 

^The SiO abundance used towards K3-50A for deriving outflow kinematics is the average abundance from the other six sources. 



Table 2. HCO"*" Observations and Derived Infall Properties 



Source 


Tinb 


FWZP'^ 


< / T*,^dv > 


<log[NHCO+]> 


<Vi,> 




Age 








(K) 


(km s~^) 


(K km s"^) 


(cm~^) 


(km s~^) 


(xlO-3 Mq yr-i) 


(xlO^ yr) 


(Mo) 




G5.89-0.39 


31.5 


140 


35.6 


13.17 






0.29 






GlO.47+0.03 


11.7 


60 


11.2 


12.67 


1.02±0.1 


8.89 


0.88 


78.0 


0.42 


GlO.6-0.4 


31.2 


50 


28.4 


13.08 


1.19±0.1 


27.34 


1.42 


389 


0.7^ 


G19.61-0.23 


8.3 


80 


14.4 


12.78 


1.49±0.6 


<8.5 


<0.71 


<60.3 


<0.3^ 


G20.08-0.14 


5.9 


55 


9.2 


12.59 


1.35±0.3 


5.51 


0.60 


33.0 


0.09 


G29.967-0.02 


16.5 


30 


9.2 


12.59 






1.19 






K3-50A 


17.0 


30 


13.2 


12.75 






<1.36 







^FuU Width at Zero Power of the HCO+ hne. 

^Infall velocity averaged over each pixel with an infall signature 
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Fig. 1. — G5. 89-0. 39 The solid contours represent integrated SiO emission towards G5.89. 
They start at 3cr (9.68 K km s"^) and increase in steps of 3a. The dashed contour is the 3a 
(5.23 K km s~^) contour of H^^CO+ emission while the spectra represent the HCO+ emission 
where the integrated intensity is above 10a (14.50 K km s~^) from -50 to 80 km s~^. 
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Fig. 2. — GlO.47+0.03 The solid contours represent integrated SiO emission towards 
G10.47. They start at 3a (1.14 K km s~^) and increase in steps of 3a. The dashed contour is 
the 3a (0.75 K km s~^) contour of H^^CO"*" emission while the spectra represent the HCO"*" 
emission where the integrated intensity is above lOcr (10.00 K km s~^) from 45 to 90 km s~^. 
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Fig. 3. — GlO.62-0.38 The solid contours represent integrated SiO emission towards G10.6. 
They start at 3cr (1.01 K km s~^) and increase in steps of 3a. The dashed contour is the 3a 
(1.13 K km s~^) contour of H^^CO+ emission while the spectra represent the HCO+ emission 
where the integrated intensity is above 10a (8.55 K km s~^) from -30 to 30 km s~^. 

10 ri ' 1 ' \ ' 1 ' 1 — ' 1 ' 1 " 




-ll'57'oo' - 

_l I I I \ 1 I I \ 1 I I I 1 

18''27"40f0 39'5 39f0 3B'5 36fO 37'5 37f0 

EA (J2000) 

Fig. 4. — G19.61-0.23 The solid contours represent integrated SiO emission towards G19.61. 
They start at 3a (1.35 K km s"^) and increase in steps of 3a. The dashed contour is the 3a 
(0.9 K km s~^) contour of H^^CO"*" emission while the spectra represent the HCO"*" emission 
where the integrated intensity is above lOcr (8.39 K km s~^) from 10 to 80 km s~^. 
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Fig. 5. — G20.08-0.14 The solid contours represent integrated SiO emission towards G20.08. 
They start at Sex (1.10 K km s"^) and increase in steps of 3a. The dashed contour is the 3a 
(0.90 K km s~^) contour of H^^CO+ emission while the spectra represent the HCO+ emission 
where the integrated intensity is above lOcr (8.23 K km s~^) from 10 to 80 km s"^. 
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Fig. 6. — G29.96-0.02 The solid contours represent integrated SiO emission towards G29.96. 
They start at 3a (1.53 K km s"^) and increase in steps of 3a. The dashed contour is the 3a 
(0.76 K km s~^) contour of H^^CO"*" emission while the spectra represent the HCO"*" emission 
where the integrated intensity is above lOcr (6.3 K km s~^) from 90 to 110 km s~^. 
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Fig. 7. — K3-50A The solid contours represent integrated SiO emission towards K3-50A. 
They start at 3a (1.15 K km s~^) and increase in steps of 3a. The dashed contour is the 3a 
(1.00 K km s~^) contour of H^'^CO^ emission while the spectra represent the HCO+ emission 
where the integrated intensity is above lOcr (5.97 K km s~^) from -50 to km s~^. 
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Fig. 8. — Spectra from G10.47 averaged over 9 pixels each. The strongest (sohd) hne 
represents the HCO+ emission towards the source center, while the dashed and dashed-dotted 
spectra represent the HCO+ emission towards the blue and red shifted regions (respectively). 

The least intense (solid) line represents the H^^CO"*" emission towards the same central region 
as the strong HCO+ line. Note that the central HCO+ emission spectrum has a double peaked 
line profile due to self absorption, while the H^^CO"*" appears single peaked and optically thin. 
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Fig. 9. — First moment map of G10.47 (gray scale). The grayscale ranges from 63 to 72 K 
km s~^. The velocity weighted integrated intensity shifts smoothly from the blue lobe to 
the central source to the red outflow lobe, suggesting that this is a single source with an 
outflow. The second moment map (black contours) shows that the highest velocity dispersion 
is centered on the location of the HIl regions. The contours start at 10% of the peak velocity 
dispersion (0.41 km s~^) and increase in 10% intervals. 



